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V 

Sand  movement  on  the  coast  by  wind  action  can  be  a  major  facto- 
in  some  areas,  yet  very  little  is  known  about  the  basic  mechanisms  of 
movement  by  wind.  With  the  growing  use  of  artificial  placement  of  sand 
material  for  beach  restoration  and  protection,  movement  of  sand  by  wind 
in  the  area  shoreward  of  the  waterline  has  become  more  and  more  a  matter 
of  concern.  The  need  for  greater  knowledge  in  this  area  has  recently 
been  emphasized  by  proposed  plans  for  the  rehabilitation  of  vast  off¬ 
shore  .and  island  areas  along  the  cu.sici.ii  seaboard.  This  report  dis¬ 
cusses  some  earlier  experimental  results  in  wind  tunnels,  and  describes 
and  compares  results  obtained  from  new  wind  tunnel  tests.  Work  by 
another  investigator,  directly  related  to  the  author's  work,  has  bei_n 
included  in  an  addendum.  This  latter  work  extends  the  investigation  to 
smaller  sand  size  range  and  suggests  that  different  relationships  ior 
threshold  shear  velocity  for  sand  movement  obtain  for  the  smaller  sard 
size  distributions.  Of  particular  interest  are  tests  on  the  influence 
of  moisture  content  on  sand  movement.  The  experimental  data  clearly 
demonstrate  that  as  the  sand  moisture  content  increases,  the  value  of 
the  threshold  shear  velocity  of  sand  movement  may  also  materially 
increase.  Quantitative  expression  of  this  effect  is  obtained.*. 

This  is  the  first  of  the  technical  memoranda  series  to  be  issued 
by  the  II,  S  Army  Coastal  Engineering  Research  Center,  the  Federal  agency 
established  to  carry'dn  the  research  function  of  the  Beach  F.tosion  Board 
which  was  abolished  on  7  November  1963.  " 

This  report  was  prepared  at  the  Wavr  Research  Laboratory  of  the 
Institute  of  Engineering  Research  of  the  University  of  California  at 
Berkeley  in  pursuance  of  contract  DA-49-055-eng-17  with  the  Beach 
Erosion  Board  which  provides  in  part  for  the  study  of  sand  movement  by 
wind.  The  author  of  this  report  Pierre-Yves  Belly,  was  a  graduate 
student  a  -  the  University  during  the  time  the  work  described  in  this 
report  was  carried  out.  During  this  time  Mr.  Belly  was  in  this  country 
on  a  French  Government  scholarship.  He  is  currently  serving  with  the 
French  Navy,  The  author  of  Addendum  II,  Abdel-Latif  Kadib,  was  also 
a  graduate  stident  at  the  University  while  this  work  was  carried  out. 


This  report  is  published  under  authority  of  Public  Lai  166,  70th 
Congress,  approved  July  ji  ,  VtVj ,  as  supplemented  by  Public  Law  172, 
88th  Congress,  approved  7  November  1963. 
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SAND  N'WEMByr  BY  NINE 
by 

Pierre-Yves  Belly 


INTRODUCTION 

Sand  movement  by  wind  action  has  already  been  treated  by  several  re¬ 
search  workers.  In  the  following  sections,  sow  facts  and  theories  related 
to  the  subject  of  sand  movement  by  wind  will  be  briefly  presented. 

Wind  velocity  above,  a  sand  surface. 

The  shear  stress,  t,  produced  at  the  sand  surface  by  wind  is  one  of 
the  most  important  factors  in  investigating  sand  movement  by  wind  action. 
When  the  aftear  stress  exceeds  a  critical  value,  the  sand  particles  start 
to  move.  When  sand  is  being  transported,  the  air  above  the  sand  surface 
behaves  as  a  heavy  and  non-bomogeneous  fluid,  so  thst  the  wind  velocity* 
distribution  is  changed,  although  the  basic  equation  remains  the  same. 

As  long  as  there  is  no  sand  movement,  the  wind  velocity  distribution 
can  be  adequately  described  by  the  general  equation, 

U  -  C  log  L. 

Zo 

in  which  U  it  the  velocity  at  height  Z  above  the  sand  surface  and  Z0  is  ■ 
reference  parameter.  The  coefficient  C,  according  to  Von  Kaxman*  s  develop¬ 
ment!^,  is  equal  to  (2.3AMJ*,  where  l  is  the  Raman  Constant,  0*  is  the 
shear  velocity  defined  as  VT/p,  and  p  is  the  density  of  air.  Taking  the 
value  of  0.40  for  I,  the  Von  Kaxman  equation  can  be  formulated  as, 

0  -  5. 75  U*  log  |- 

Concerning  the  roughness  factor,  ZQ,  Zingg^  proposes  the  equation, 

Z0  -  0.081  log  1J~ 

with  Z0  and  the  sand  grain  diameter,  d,  expressed  in  millimett  a,  This 
equation  contains  both  the  results  of  Bagnold'2^  (Zo  *  d/30)  for  aaall 
train  sizes,  and  that  of  White^  (Z0  *  d/9)  for  large  grain  sizes. 

Once  the  wind  velocity  is  great  enough  to  move  sand  particles,  the 
wind  velocity  distribution  is  altered  by  the  sand  movement.  Plotted  on 
semi-log  paper,  the  velocity  distributions  remain  straight  lines,  but  as 
shown  by  Bagnold,  they  ill  seem  to  meet  -t  a  certain  poiut,  which  he  called 
"a  focus".  The  height  of  the  focus,  Z’,  appears  to  be  associated  with  the 
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height  of  the  ripples  which  fora  on  the  surface,  in  t  way  somewhat  analo¬ 
gous  to  that  in  which  Z0  is  associated  with  the  dimensions  of  the  grains. 
For  a  sand  of  average  grain  size  0.25  nm;  Bagnold  found  tint  height  of  the 
focus,  Z*t  to  be  about  3  an,  a td  the  corresponding  velocity  to  be  about 
2.5  it/ sec.  A  more  thorough  study  wade  by  Zingg^4*  allows  one  to  predict 
the  focus  by  means  of  the  formulae, 

Z’  v  10  d  millimeters 

U’  >  20  d  miles/hour 

where  the  grain  diameter,  d,  is  expressed  in  millimeters.  Using  the  com¬ 
ponents  of  the  focus,  Z*  and  U* ,  the  wind  velocity  distribution  can  be 
expressed  by, 

U  *  C  log  ♦  u* . 

Bagnold  assianed  e  coefficient  C  of  5.75  U*,  which  corresponds  to  the  value 
o'  .'.*•>  for  the  Ktrman  Constant.  But  experiments  by  Zinggf^J  yielded  the 
equation, 

U  -  6.13  U*  log  |,  *  U* 

which  indicates  a  value  of  0.375  for  the  Karman  Constant. 

Sand  movement  by  wind. 

When  the  wind  above  a  sand  surface  is  great  enough,  the  particles 
start  moving.  The  wind  velocity  profile  and  the  shear  velocity  age  the 
primary  motivating  factors  in  initiating  and  sustaining  sand  movement. 

The  initiation  of  sand  movement  has  been  theoretically  investigated  by 
Bagnold^).  He  obtained  for  the  threshold  value  of  the  shear  velocity. 


where  d  is  the  grain  diameter,  g  is  the  acceleration  of  gravity,  a  and  p 
are  specific  weights  of  sand  and  air,  respectively.  Bagnold  found  that  the 
coefficient  A  is  nearly  constant  for  a  sand  diameter  of  0.25  mt  and  for  all 
sands  of  larger  grain  size.  A  approximates  0.1.  Experiments  by  Zingg 
corroborate  this  result.  For  very  wall  grains  (when  the  Reynolds  number 
u*  d/v  is  less  than  the  critical  value  3.5,  i.e.  when  the  surface  becomes 
"smooth")  the  value  of  the  coefficient  A  is  no  ljnger  constant.  Figure  1 
shows  the  variation  of  the  threshold  velocity  with  grain  size  as  found  by 
Bagnold. 


(0-: 


U»t  in  cm/sec 


30  \ 


Groin  Diameter  in  mm 


,01.04  il  .2  .4  ,6  .8  I  1.5  2  2  5 


F'GURE  i-  INITIATION  OF  SAND  MOVEMENT 


The  theories  on  sand  movement  can  be  classified  Into  two  groups.  One 
is  based  on  the  investigation  of  the  vertical  distribution  of  the  sand  move¬ 
ment  above  the  bed,  and  the  other  is  based  on  the  assiasption  that  the  sand 
particles  move  downstream  with  bouncing  motions  near  the  sand  surface. 
Representative  of  the  former  concept  are  theories  of  Kawamura  and  of  Ishi- 
hara  and  Iwagaki  (see  reference  1);  experimental  results  correspond  well 
with  the  above  theories,  and  the  total  rate  of  transport  could  be  obtained 
by  integration  with  respect  to  the  height,  but  the  expressions  are  too 
complicated  for  practical  use.  The  Bagnold  and  Kswamura  formulae  are 
obtained  by  using  the  latter  theories,  and  are  expressed  below. 

,  (2) 

Bagnold  formula 


The  rate  of  sand  movement  per  unit  width  and  unit  time,  q,  is  given  by. 


q " c  Jr 


where  D  is  the  grain  diameter  of  a  standard  0.25  asa  sand,  d  is  the  grain 
diameter  of  the  sand  in  question,  p  is  the  specific  weight  of  the  air 
(p/g  *  1.25  10-t>  e.g.s.),  04  is  the  shear  velocity  and  C  has  the  following 
values: 

1.5  for  a  nearly  uniform  sand 

1.8  for  a  naturally  grided  sand 

2.8  for  a  sand  with  a  very  wide  range  of  grain  sizes. 


1- 


tawaaura  formula^ 3  ^ 

The  rate  of  sand  movement,  q,  is  given  by, 

q  =  k  |  (U*  -  Un>  (U*  ♦  U*t>2 

where  p  i*  the  specific  weight  of  air,  U*  is  the  shear  velocity,  is  the 
threshold  shear  velocity,  and  K  is  a  constant  which  should  he  determined  by 
experiment.  For  a  sand  of  average  grain  sire  0.25  nm,  Kewamura  obtained 
k  «  2.78  in  a  wind  tunnel. 


The  basic  ideas  of  the  above  formulae  are  almost  the  same,  and  excepr 
for  light  winds,  both  relationships  give  approximately  the  same  results 
if  a  suitable  constant  is  chosen.  On  the  contrary,  as  shown  by  Figure  2, 
experimental  results  obtained  in  wind  tunnels  by  Bagno  Id  and  Kawatrura 
differ  widely,  although  the  sand  diameter  was  almost  the  same,  Experimental 
results  obtained  in  wind  tunnels  ’oy  Zingg(4)  and  HorUtswi.i  1)  are  also 
plotted  on  Figure  2.  From  his  results  Zingg  modified  the  Bagno  Id  formula 
thus. 


with  C  -  0.83 


0  . 


d  S5/4 


? 


g 


f  8) 

In  addition  to  these  theoretical  formulae,  O'Brien  and  Rindlaub 
proposed  tne  following  formula  from  data  derived  in  the  field: 


G  ■  0.036  U53  <for  U5*20  ft/see) 


where  G  ii  the  rate  of  transport  in  pounds  per  dey  per  foot  width,  and  U 
is  the  wind  velocity  5  ft,  above  the  sand  surface  in  ft/sec. 


Confirmation  of  these  formulae  by  field  results  is  not  particularly 
good,  but  since  there  is  considerable  scatter  in  the  experimental  data, 
these  formulae  are  still  useful  in  the  description  of  a  particular  con¬ 
dition  when  a  suitable  constant  is  chosen. 

Flying  distance  of  sand  particles. 

(2  ) 

Mathematical  approaches  to  this  problem  have  been  made  by  Bagnold 
and  Ford("\  Photographic  observation  of  the  sand  path  confirmed  theo¬ 
retical  results  in  both  cases.  From  the  sand  distribution  in  a  horizontal 
trap,  KawamuraO)  and  Horibawaf*)  derived  the  average  flying  distance  of 
the  sand  particles.  Experimental  values  are  plotted  in  Figure  3. 

Sand  traps. 

Experiments  by  Horikawa  and  Shen^  demonstrated  that  none  of  the 
available  sand  traps  used  by  previous  investigators  gave  entirely  satis¬ 
factory  results.  They  showed  that  the  efficiency  of  a  horizontal  trap  can 
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FIGURE  5.  WIND  TUNNEL 


a 


be  aide  relatively  high,  simply  by  making  the  trap  reasonably  long.  As 
for  the  vertical  type  of  sanu  trap,  they  developed  one  in  which  the  dis¬ 
turbance  of  the  flow  is  minimum,  the  efficiency  approaches  100%,  and  it  can 
be  used  for  experiments  in  wind  tunnels  as  well  as  for  field  experiments. 
This  vertical  trap  is  described  later. 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Experiments  were  conducted  in  a  wind  tunnel  located  in  Building  276 
at  the  Richmond  ?ield  Station  of  the  University  of  California.  This  tunnel, 

4  ft.  wide,  2.5  ft.  high,  and  100  ft.  lone,  was  constructed  cf  olywood;  the 
lower  part  of  one  side  was  made  of  glass  for  observation  (Figures  4  and  *>}. 
The  wind  was  generated  by  a  fan  at  the  exit  end.  The  mean  velocity  was 
varied  from  24  to  40  ft/sec  by  a  rheostat  controlling  the  fan  speed. 

Wind  velocities  were  measured  with  a  standard  Prandtl  type  pitot  tube 
which  was  attached  to  a  point  gage  and  introduced  into  the  air  stream 
through  the  top  of  the  flume.  The  pitot  tube  was  connected  to  a  Magnehelic 
gage  having  a  range  of  one-half  inch  of  water  and  graduated  into  divisions 
of  0.02  inch.  The  Magnehelic  gage  was  chosen  instead  of  an  Ellison  draft 
gage,  because  its  response  to  pressure  changes  is  more  rapid. 

The  mean  diameter  of  the  sand  used  in  this  experiment  was  0,44  mm,  a  f 
shown  by  the  mechanical  analysis  curve  in  Figures  6a  and  6b.  The  sand  was 
spread  over  a  length  of  62  ft.  of  the  fluate,  with  a  thickness  of  2  inches. 

A  hopper  was  placed  at  the  entrance  of  the  flume,  for  use  during  long  runs, 
and  the  feed-in  adjusted  to  the  rate  of  transport. 

The  vertical  trap  developed  by  Horikawa  and  Shent1^  was  placed  in  the 
center  of  the  flume  7  ft.  before  the  end  of  the  sand  bed.  This  trap  had 
a  width  of  3/8  inch  and  a  height  of  1  ft.  (Figure  9).  A  horizontal  trap 
8  ft.  long,  consisting  of  18  compartments  was  permanently  fixed  at  the  end 
of  the  sand  bed.  In  order  to  avoid  the  side  wall  effect  the  amount  of  sand 
retained  in  the  trap  was  measured  only  in  the  center  part  of  the  flume  over 
a  width  of  2  ft. 

The  de^'.zed  wind  velocity  was  obtained  by  adjusting  the  rheostat  of  the 
speed  control  on  the  fan.  Each  run  was  allowed  to  continue  for  a  period  of 

5  to  30  minutes.  The  ripples  on  the  bed  and  the  scour  around  the  vertical 
trap  were  observed.  At  the  end  of  each  run  the  horizontal  trap  was  cleaned 
out  with  a  vacuum  cleaner  (Figure  4). 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 
Velocity  distribution  on  sand  surface 

In  order  to  investigate  the  side  wall  effect  on  wind  velocity,  the  wind 
velocity  distributions  in  a  cross  section  were  measured.  The  transverse 
profiles,  11  ft.  before  the  end  of  the  sand  bed  are  shewn  in  Figure  10  for 
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different  fan  currents.  These  transverse  profiles  show  tne  wind  to  be 
practically  unifora  across  the  channel  except  in  the  close  proximity  of 
the  walls. 

Vertical  wind  profiles  were  measured  at  the  same  place  (11  ft.  before 
the  end  of  the  sand  bed)  at  the  center  of  the  flume.  Tne  velocity  dis¬ 
tributions  obtained  for  different  fan  currents  are  shewn  in  Figure  11,  and 
plotted  on  seal-log  scale  in  Figures  12  and  13. 

For  wind  velocities  less  than  the  critical  value  required  to  initiat- 
sand  movement,  the  relationship  between  wind  velocity  and  height  above  the 
sand  surface  obeys  the  logarithmic  law  (Figure  12). 

For  wind  velocities  larger  than  the  critical  value,  the  relationship 
also  obeys  the  logarithmic  law  above  the  focal  point  (Figure  13).  The 
focal  point  located  at, 

Z*  ■  0.0144  ft. 

U*  1 3  ft/ sec 

seems  to  agree  with  Zingg* s  estimate  of, 

Z'  =  10  d  am 

U*  »  20  d  miles  per  hour 

where  d  is  the  wan  diameter  of  the  sand  in  millimeters,  or: 

Z'  *  0.0135  ft. 

IT  -  13  ft/sec 

The  shear  velocity  U*  can  be  determined  by  the  slope  of  the  velocity 
distributions  in  Figure  13.  The  values  of  for  different  wind  velocities 
(measured  at  Z  *  1.0  ft.  above  the  sand  bed),  calculated  from  Zingg' s 
formula , 

U  =  6,13  U*  log  |,  *  O', 


are  given  in  Table  1. 
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Table  1 


Wind  velocity  Shear  velocity 

{at  Z  «  1.0  ft.)  in  cm/sec 

in  ft/sec 


25.0 

25.7 
26.0 
27.0 
28.2 
30.0 
31.0 

32.8 
34.5 
39.0 


39.0 

37.5 

42.2 
43.8 

46.3 

47.6 
51.0 

55.5 

58.5 
70.0 


The  relationship  between  U  and  U4  is  approximately  linear  (Figure  14). 
The  shear  velocity.  { }*,  for  other  wind  velocities  was  determined  by  using 
this  graph. 

The  value  of  the  threshold  velocity  is  radically  changed  by  the  pres¬ 
ence  or  the  absence  of  a  sand  feed-in.  Without  sand  feed-in.  0^  a  40 
cm/sec ,  but  with  sand  feed-in  the  threshold  velocity  is  greatly  lowered: 

U*t  m  30  cm/sec.  The  latter  value  is  very  close  to  the  value  calculated 
with  the  Bagnold's  formula  (0*t  «  34  cm/s ec). 

Rate  of  sand  transport 


The  amount  of  sand  caught  by  the  horizontal  trap  was  measured  for 
velocities  varying  from  the  threshold  value  to  37  ft/sec.  The  sand  feeding 
which  was  located  at  the  upstream  end  of  the  sand  bed  is  a  very  important 
factor  in  sand  movement  for  lower  wind  velocities,  The  sand  feeding 
greatly  lowers  the  threshold  velocity  and  at  the  same  time  changes  the 
amount  of  sand  transported  for  lower  velocities.  Table  2a  gives  the  data 
obtained  without  sand  feeding.  Table  2b  givea  the  results  obtained  when 
the  sand  feeding  was  established  with  a  discharge  approximating  the  rate 
of  sand  transport.  These  results  are  plotted  in  Figure  15  and  in  Figure 
16  to  show  the  comparison  with  those  of  Bagnold,  Kawamura  and  IJorikawa. 

The  apparent  reversal  of  the  curve  obtained  without  sand  fceuiiig  is 
perhaps  due  to  the  fact  that  the  sand  used  in  this  study  has  a  wide  range 
of  grain  sizes  (0.2  to  0.7)  mm).  At,  or  near  the  threshold  it  is  possible 
that  the  action  of  the  smaller  grains  was  impeded  by  the  larger,  thus 
modifying  the  over-all  values  for  the  threshold  and  rate  of  transport. 

More  precisely,  near  the  threshold  t>e  st  d  grains  move  mainly  by  saltation. 
Since  the  surface  layer  remains  practically  immobile  (no  surface  creep). 
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the  smaller  grains  are  hidden  by  the  larger  ones  and  as  j  result  the  sand 
behaves  as  if  it  had  a  much  larger  mean  diameter.  According  to  the  Bagnold 
formulae  for  the  rate  of  transpot  t  arm  the  th.eshc  ld  value  of  the  shear 
velocity,  the  curves  for  too  different  mean  riameters  arc  as  sketched  below; 


Therefore  the  sand  which  initially  follows  curve  (1),  gradual!,  changes 
its  effective  mean  diameter  and  begins  to  follow  curve  (2).  This  phe 
nomena ,  which  is  related  to  the  state  of  the  surface  layer,  disappears 
when  this  surface  layer  is  artificially  set  in  motion  by  the  sand  feed-in. 
Therefore  no  anomaly  is  noted  in  the  resulting  curve. 

The  experimental  values  for  the  maximum  rate  of  transport  (i.e.  with 
sand  feed-in),  q,  can  be  compared  to  the  values  predicted  from  the  Bagnold 
and  Kawamura  formulae. 

For  the  average  grain  sire,  d  *  0.44  *m,  the  Bagnold  formula  gives  in 
e.g.s.  units, 

I o'  44  -6  3 

q  =  c  Voi  X  1.25  x  10  x  U* 

Taking  C  -  2.5  (Bagnold  proposes  1.8  for  normally  graded  sand  and  2.3  for 
a  aand  of  a  very  wide  range  of  grain  size),  this  formula  is  plotted  in 
Figure  17.  Except  for  wind  velocities  approaching  the  threshold,  which  in 
any  case  cannot  be  described  by  Bagnold '$  formula,  the  agreement  with 
experimental  results  is  very  good. 

The  Kawamura  formula  is  in  e.g.s.  units, 

q  *  k  x  1.23  x  10"6  x  (U#  -Utt)  (U#  +  U*t>2 
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With  sand  feed- in  we 

found  U.t  to  be  about  30  ca/atc. 

| 

Putting  this  value 

in  the  Kawaaura  formula  and 

using  X  *  3.1f  *hia  fnr*u2*  describes  very  well 

the  rate  of  sand  transport 

for  the  whole  i*nt:e  of  t! 

ovities  (Figure  18). 

Thus,  by  givinp  to  the  constants  adequate  value*, 

t>*  formulae  of 

Bagnold  and  Kawamura 

agree 

very  well  with  the  result;! 

obtained. 

Tabic  ?.m 

Kind  velocity 

Shear  velocity 

P.ate  of  transport 

(at  Z  *  1.0  ft. ) 

in  cm/sec 

in  gr/cm-sec 

in  ft/sec 

25.0 

39 

none 

26.0 

41 

0,143 

26.0 

41 

0.012 

26.2 

42 

0.066 

26.8 

43 

0.137 

27.0 

43 

0.096 

37  0 

43 

0.303  j 

27.2 

44 

0.187 

27.8 

45 

0.313  i 

27.9 

45 

0.292 

28.3 

46 

0.386 

28.5 

46 

0.382 

30.5 

50 

0.505 

30.3 

51 

0.545 

31.0 

52 

0.580 

32.7 

54 

0.700 

33.0 

55 

0.780 

34.8 

58 

0.910 

37.0 

64 

1.18 

Table  2b 

Wind  velocity 

Shear  velocity 

Rate  of  transport 

(at  Z  »  1.0  ft.) 

in  cm/sec 

in  gr/cm-sec 

in  ft/sec 

20.0 

30 

0.012 

23.0 

35 

0,105 

24.8 

38 

0.182 

25.0 

39 

0.220 

26.0 

<1 

0.232 

28.4 

46 

0.380 

34.0 

50 

0.506 

17 

.US  Jr 


S». 

^■1 


FIGURE  14  -  RELATIONSHIP  BETWEEN  THE  MAIN 
VELOCITY  U,  AND  THE  SHEAR  VELOCITY  U*. 
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Efficiency  of  the  vertical  sand  trap 

The  efficiency  of  the  vertical  sand  trap  was  tested  for  various 
velocities  in  the  course  of  experiments  on  the  rate  of  transport.  The 
horizontal  trap  was  long  enough  to  catch  practically  all  of  the  sand 
transported  by  the  wind  and  served  as  a  reference  for  the  vertical  trap. 
The  efficiency  of  the  vertical  sand  trap  is  defined  as. 


„ _  amount  of  sand  caught  in  vertical  trap  ... 

11%  ■  '  '  .  -u-  ""  ■  ■  . - ! - ;  'V'".  r  "  x  100 

amount  of  sand  caught  in  horizontal  trap 

Table  3  and  figure  19  give  the  efficiency  for  different  wind  velocities. 

Table  3 


Wind  velocity 

Vert,  trap 

Her.  trap 

Efficiency 

in  ft/sec 

q  in  gr/cm-sec 

q  in  gr/cm-sec 

% 

25.8 

0.084 

0.143 

60 

26.0 

0.004 

0.116 

38 

27.0 

0.26 

0.30 

87 

27.0 

0.022 

0.096 

23 

27.8 

0.267 

0.313 

85 

28.3 

0.340 

0.386 

88 

28.5 

0.314 

0.382 

83 

30.5 

0.48 

0.505 

93 

31.0 

0.53 

0.58 

92 

31.5 

0,51 

0.51 

100 

31.5 

0.48 

0.50 

105 

32.0 

0,65 

0.69 

107 

32.0 

0,58 

0.61 

106 

35.0 

0.90 

0,98 

110 

35.0 

1.01 

1.33 

112 

35.0 

0.94 

1.02 

108 

35.5 

0.98 

1.13 

118 

i 


An  efficiency  higher  than  100%  probably  is  due  to  the  small  amount  of 
unmeasured  sand  which  fell  beyond  the  horizontal  trap,  and  also  to  the 
possibility  of  secondary  currents  in  the  proximity  of  the  mouth  of  the 
vertical  trap.  In  any  case  the  vertical  trap  has  a  sufficiently  good  effi¬ 
ciency  for  velocities  between  30  and  33  ft/sec  to  avoid  the  necessity  of 
corrections  in  later  experiments. 


Almost  immediately  after  the  beginning  of  a  run,  a  scour  takes  place 
around  and  below  the  vertical  trap  as  shown  on  Figure  8,  This  scour  seems 
to  remain  steady  and  therefore  does  not  influence  the  measurements.  For 
runs  of  long  duration,  however,  the  platform  becomes  undermined,  and  this 


21 


is  probably  a  esuse  of  error,  for  the  surface  creep  does  rot  thereafter 
enter  the  mouth  of  the  trap,  as  sketched  below: 


This  phenomenon  which  occurs  for  a  run  duration  of  about  one  hour  was 
avoided  as  much  as  possible  by  using  run  times  of  5  to  15  minutes  (except 
in  the  last  part  of  the  experiment  where  runs  had  durations  of  30  to  45 
minutes). 

For  higher  velocities,  the  grain  sire  distribution  of  the  sand  caught 
in  the  vertical  trap  is  very  close  to  the  grain  sire  distribution  of  the 
bed  (Figure  20).  The  relative  absence  of  bigger  grains  in  the  sand  caught 
by  the  vertical  trap  is  probably  caused  by  the  platform  of  the  trap  which 
can  sometimes  be  an  obstacle  to  surface  creep.  For  velocities  approaching 
the  threshold  value,  the  grain  sire  distribution  in  the  vertical  trap  shows 
a  distinct  lack  of  the  larger  grains  (Figure  21).  This  fact  cannot  be 
entirely  attributed  to  the  inefficiency  of  the  vertical  trap  and  probably 
indicates  the  manner  in  which  the  sand  is  moving  neat  threshold  (the  large 
grains  not  taking  part  in  the  general  movement). 

Ripples  on  sand  surface 


The  ripples  produced  on  the  sand  surface  were  observed  during  the 
different  parts  of  the  experiment.  They  appear  on  a  flat  sand  surface  as 
soon  as  there  is  some  sand  movement  and  they  disappear  at  very  high  ve¬ 
locities  (about  36  ft/sec)..  The  wave  lengths  of  the  ripples  were  measured 
for  different  wind  velocities  but,  as  shown  in  Figure,  22,  there  is  no  clear 
correspondence  between  wave  length  and  wind  velocity.  The  average  wave 
length  is: 

\  •  3  inches. 


ZZ 


Percentage  Percentage 


FIGURE  20 -GRAIN  SIZE  DISTRIBUTION  FOR  U>27  ft/sec 


FIGURE  21  -GRAIN  SIZE  DISTRIBUTION  FOR  U=25  ft/sec 
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X  in  inches 


No 

h  Sand 
Movement 


Wave  length  of  the  ripples 


<*>0® 


o 


o 


No 

Ripples 


U  in  ft  /sec 


125  25 30  35  40 

FIGURE  22 -WAVE  LENGTH  OF  THE  RIPPLES 


FIGURE  23-  RIPPLES  ON  SAND  SURFACE 
U  1  31  ft/sec 
X  *  27  inches 

H  8930 

25 
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Fig.  24 

The  x-axis  is  tiken  in  the  direction  of  the  wind.  If  30  is  tb?  Mount  of 
sand  felling  on  a  unit  area  of  sand  surface  during  unit  tiae,  the  Mount 
of  sand  "jumping"  frosi  the  surface  of  unit  width  and  of  length  dx  is: 
vjy.dx.  If  we  acsuae  now  that  all  the  particles  he.ve  the  mac  flying  dis¬ 
tance,  L,  the  sand  particles  which  pass  over  the  section  (O)  are  juaping 
froa  the  sand  surface  between  the  sections  (C j  and  (A),  with  OA  a  L. 
Therefore  the  aaount  of  sand  passing  through  the  cross  section  (0)  per 
unit  tine  is, 

«  -  £  G0.d*  -  G0.L 

Thus  the  sversge  flying  distance  is. 


i  ■  3" 

Go 

q,  the  rate  of  transport  has  been  already  obtained  by  the  experiaent 

Gq,  which  nay  be  considered  as  the  aaount  of  sand  falling  per  unit 
width  and  unit  tiae,  in  a  trap  of  infinitely  sasll  length  placed 
iaaediately  after  the  end  of  the  sand  bed,  can  be  obtained  by 
extrapolating  the  curve  of  sand  distribution  in  the  horizontal 
trap  for  x  a  0, 

As  the  first  trap  has  a  length  of  0,5  foot,  the  percentage  of  sand 
fro#  the  surface  creep  is  very  snail,  and  therefore  it  is  unnecessary  to 
correct  the  curve  of  sand  distribution,  so  that  the  aaount  of  aand  trans¬ 
ported  by  saltation  can  be  blended  with  the  total  rate  of  transport,  q. 

Froa  Figures  25,  26  and  27  we  have: 

For  U  •  28  ft/sec,  (^  ■  j8  lb/ft2-houi ,  and  L  •  4--  «  1.30  ft. 


26 
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20  25  30  35 

Wind  Velocity 

FIGURE  28 -AMOUNT  OF  SAND  FALLING  ON  UNIT 
AREA  OF  SAND  BED  PER  UNIT  TIME,  G0 


W-nd  Velocity 

FIGURE  29  -AVERAGE  FLYING  DISTANCE 


For  U  *  31  ft/sec ,  G  «  90  lb/ft2-hour,  and,  L  *  »  1.35  ft. 

For  0  ■  35  ft/aec,  G0  ■  140  lb/ft2-hour,  and  L  ^  ■  l.J  ft. 

Figures  23  and  29  give  Go  and  the  average  flying  distance,  L,  as  a  function 
of  wind  velocity. 

For  the  sane  values  of  the  velocity  but  for  a  sand  of  0.25  mm,  Bagno Id 
found  respectively  2,8,  3  and  5  inches,  for  the  average  flying  distance. 

The  average  flying  distance  increases  with  the  grain  disaster  of  the  sand, 

but  more  study  shou be  done  to  determine  if  values  as  large  as  found  ir, 

these  experiments  sre  ressonsble.  Also,  Bagnold  found  s  recarkable  agree¬ 
ment  between  the  average  flying  distance  and  the  wave  length  of  the  ripples 
produced  on  the  sand  surface.  As  the  wave  length  in  the  present  tests  was 
a  constant  value  of  3  inches,  the  Bagnold  relationship  was  not  verified. 

Frequency  distribution  function  of  flying  length 

First,  consider  the  condition  as  shown  in  Figure  30 


The  x-axis  is  taken  in  the  direction  of  the  wind  and  the  region  of  x<0  is 
covered  with  sand,  from  where  the  sand  particles  are  flying  and  dropping 
into  the  sand  trao  set  at  x>  0,  The  amount  of  sand  falling  in  the  unit 
arcs  ir.  x*  0  per  unit  time,  P  (x),  it  the  summation  of  sand  particles 
jumping  from  the  sand  surface  in  the  region  x<  0.  Therefore  if  g(L>  is 
the  frequency  distribution  function  of  the  flying  length.  Fix!  is  expressed 
by  the  following  integration. 


F(x>  m  G0 


-?  l.d? 


where,  as  seen  before.  Go  is  the  amount  of  sand  falling  on  the  unit  ares 
of  sand  surface  during  unit  time. 
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If  t  •  x  -  5 ,  Fix)  can  b«  rewritten  such  as. 
Fix)  »  -  Ge  J  git)  dc 


or. 


Fix; 


-  0o  |  tit)  dt 


By  definition, 

-  F  git)  dt  m  l 

Therefore, 

Fix)  -  Ge  [1  -  j  g(t)  dt]  . 


Then,  the  distribution  function  of 
Fix),  with  respect  to  x, 

gix)  «  - 


flying  length  is  obtained  by  differentiating 

1_  dFix) 

G0  *  dx 


The  functions: 

Fix)  •  G0  e  '  0,390  x  ,  for  U  ■  2d  xt/see 

Fix)  «  Go  e  '  0-375  x  ,  for  V  •  31  ft/aec 

and  P(x)  =  G0  e  *  0-333  X  ,  for  U  ■  35  ft/sec 

with  x  in  feet,  have  been  found  to  fit  reasonably  the  horizontal  distribu¬ 

tion  of  sand  drift  given  in  Figures  25,  26  and  2?.  By  differentiation  with 
respect  to  x,  we  find, 

gix)  «  -  0.390  e  "  0-39w  x,  for  U  -  28  ft/sec 
-  0  375  x 

g(x)  *  -  0*375  e  *  #  for  U  ■  31  ft/sec 

gix)  «  -  0.333  e  ‘  0-333  x,  for  U  -  3S  ft/aec 
These  distributions  are  plotted  in  Figure  31. 


Variation  of  the  flying  distance  with  grain  size 

This  study  can  be  Bade  by  knowing  the  grain  size  distribution  in  the 
different  cawpsr taents  of  the  horizontal  sand  trap.  This  analysis  has 
been  done  for  the  wind  velocity  of  31  ft/avc. 
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Table  4  gives  the  weight  of  sand  (expressed  in  pounds  per  hour  and 
per  foot  width)  for  each  compartment  of  the  horizontal  trap,  and  for  each 
range  of  grain  sizes.  The  results  are  plottrj  in  Figure  32. 


Extrapolating  the  curves  for  x  *  O,  the  amount  of  sand  falling  on  a 
unit  area  of  sand  surface  during  unit  time,  Go,  can  be  netermined,  and  the 
average  flying  distance,  L,  calculated  (I.  a  q/G0)  for  the  different  ranges 
of  grain  sizes  (Table  5). 


As  a  final  result.  Figure  33  shots  the  variation  of  the  flying  dis¬ 
tance,  L,  with  grain  size. 


Table  4 


T 


U  a  31,3  ft/sec 

Distribution  for  different  grain  size  A 


Dist. 

200- 

300 1* 

300-350^ 

350 

-450H 

450 

-500H 

500- 

-600d 

600- 

-lOOOd  L 

in  ft. 

% 

q 

% 

q 

% 

q 

% 

q 

% 

q 

X 

q 

q 

0.5 

8.5 

3 

14.5 

5.1 

19 

6.6 

30 

10  5 

17 

6 

10 

3.5 

35 

1.0 

3.8 

5.9 

4.2 

4 

1 

0 

24 

1.5 

28 

4.2 

37 

5.6 

21 

3.1 

11 

1.65 

2 

0.3 

15 

2.0 

3.5 

3.0 

2.3 

0.8 

0.1 

10 

2.5 

30 

2.2 

26 

2.0 

31 

1.8 

11 

0.8 

1 

0.075 

Eli 

| 

7.5 

3.0 

1.8 

1.2 

1.2 

0.5 

0.05 

1 

- 

5 

3.5 

36 

1.45 

28 

1.1 

23 

0.92 

to 

0.35 

1 

0.035 

K 

I 

4 

4.0 

1.3 

0.7 

0.7 

0.1 

0 

2,5 

5.0 

62 

1.25 

25 

0.5 

11 

0.22 

15 

0.03 

0 

0 

2 

7.0 

IP 

EB 

8 

0.16 

0. 

02  0 

0 

0 

2 

8 

Eli 

D 

0.06 

0. 

7  0.007 

0 

0 

| 

1 

Total 

q 

24.8 

24.8 

- - 

21.7 

19.0 

7.4 

4.0 

Note:  The  amount  of  sand  "q"  are  in  lb/ft/hr 


The  values  of  q  which  could  not  be  found  in  the  experimental 
results  have  been  determined  from  the  curves  (Figure  32), 


FIGURE  31  -  FREQUENCY  DISTRIBUTION  FUNCTIONS  OF  FLYING  LENGTH 


FIGURE  3£  -  SAND  DISTRIBUTION  IN  THE  HORIZONTAL  TRAP  FOR  DIFFERENT  GRAIN  SIZES 


amount  of  *ond 

t  collected  every 
2  minutes,  in  grams 
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FIGURE  34 

N.  B.  The  overoge  of  oil  the  voiues  obtained  on  the  amount  of  sand  collected 
for  o  particular  velocity,  hos  been  taken  as  a  reference  In  plotting 
the  discrepancies.  This  hos  beer  done  in  order  to  compare  the 
magnitude  of  the  discrepancies  to  the  absolute  values  of  the  amount 
of  sand  collected,  but  has  no  quantitative  volue. 


34 


/ 


Table  3 

Flying  length: 


Grain  size  in  p 

200-300 

300-350 

350-450 

450-500 

500-600 

600-1000 

q  (lb/ft/hf) 

24.8 

24.9 

21.7 

19.0 

7.4 

IfinffRVHMi 

5 

9 

16 

28 

32 

L  (ft) 

5 

2.8 

1.35 

0.68 

0.23 

0.1 

(G  *  2  x  extrapolated  value  (obtained  for  a  1/2  ft.  long  trap)) 


Response  tine  of  the  sand-bed  to  a  change  of  wind  velocity 

In  order  to  investigate  the  response  time  of  the  bed  to  a  change  of 
wind  velocity,  the  wind  was  first  allowed  to  blow  over  the  sand  surface 
for  a  relatively  long  time  (sufficiently  long  to  observe  a  duplication  of 
the  results  on  the  amount  of  sand  transported,  both  in  the  vertical  and 
horizontal  traps).  The  wind  velocity  then  was  suddenly  changed  to  a  higher 
value,  and  the  sand  collected  in  the  vertical  trap  was  weighed  every  two 
minutes,  until  a  duplication  of  results  was  observed,  the  wind  velocity 
then  was  again  adjusted  to  the  previous  value,  while  the  same  measurements 
were  made  at  the  vertical  sand  trap.  After  a  sufficiently  long  time,  the 
wind  velocity  was  again  adjusted  to  the  higher  value,  and  the  same  process 
repeated.  The  two  particular  values  of  the  velocity  were: 

a  w  31,5  ft/sec 

U  a  35  ft/sec 

(measured  at  Z  *  1.0  ft,  above  the  sand  surface).  The  results  from  these 
tests  are  in  Tables  6  and  7.  There  is  a  considerable  scatter  especially 
for  the  higher  velocity.  This  dispersion  probably  is  due  to  some  extent 
to  the  inaccuracy  of  the  wind  velocity  readings,  these  velocities  being 
slightly  different  in  corresponding  runs.  But  as  the  important  fact  is 
the  rate  of  transport  with  respect  to  time,  we  can  eliminate  the  part  of 
the  dispersion  due  to  differences  in  the  main  wind  velocity,  by  considering 
the  discrepancy  between  the  measurements  made  within  the  first  few  minutes 
of  each  run,  and  the  average  of  the  last  measurements,  when  equilibrium  is 
reached.  Figure  34  shows  the  rate  of  sand  transport  as  a  function  of  time 
with  this  correction. 

The  dispersion  is  greater  at  the  beginning  than  at  the  end  of  a  run, 
but  except  for  a  slight  increase  during  the  first  4  minutes  in  most  of  the 
runs,  no  clear-cut  tendency  can  be  observed. 


4= 


35 


Table  6 


The  first  velocity  was  31  ft/aee,  and  the  second  was  35  ft/sec.  The 
figure*  represent  the  weight  of  sand  collected  in  the  vertical  trap, 
every  two  minutes,  in  grans. 


Tine  in  minutes 


Run  No. 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

1 

86.5 

94 

90.5 

100 

100.5 

94 

99 

3 

86 

94 

97 

96 

97.5 

99 

102 

5 

75 

79 

77.5 

78.5 

SO. 5 

84 

84 

73 

80 

71 

7 

77 

70.5 

67 

69 

70 

69.5 

9 

70.5 

74 

69 

73 

11 

102 

96 

97 

94 

91 

13 

82 

79 

90 

80 

90 

Table  7 

The  previous  velocity  was  35  ft/sec,  the  new  one  is  31.5  ft/sec 


Tine  in  minutes 


Run  No. 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

2 

47 

49 

51 

54 

r  , 

51.5 

51 

4 

46.5 

48.5 

48 

44 

50 

50.5 

51.5 

6 

44.5 

42.5 

41 

36 

38 

35 

37 

43 

36 

36 

8 

50 

32 

43 

44 

40 

45 

41 

40 

38 

41 

35 

35 

43 

38 

40 

12 

41 

49 

46 

46 

14 

28 

29 

32 

30 

31 

Additional  runs  were  made  in  order  to  detect  a  more  subtle  develop¬ 
ment.  In  those  runs,  the  sand  was  collected  at  the  vertical  trap  every 
30  seconds: 

Table  8 

First  velocity:  31  ft/sec 
Second  velocity:  35  ft/sec 

Tine  in  seconds 


Run  No. 

30 

60 

90 

120 

15 

26 

25 

26 

25 

17 

19 

19 

24 

24 

36 


Table  9 


First  velocity:  35  ft/sec 

Second  velocity:  31  ft/sec 


Tine  in  seconds 


Run  No. 

30 

L  f\ 

uu 

90 

120 

16 

10 

12 

9 

10 

IS 

6 

7 

6 

7 

Again  it  was  possible  to  detect  immediately  a  noticeable  change  in 
the  rate  of  transport,  but  no  further  development.  In  conclusion  it  can 
be  said  that  the  sand-bed  adjusts  itself  almost  immediately  to  a  new  wind 
velocity. 


SUMMARY  AND  CONCLUSIONS 


1.  Although  the  grain  size  of  the  sand  used  was  very  different  from 
that  used  by  Bagnold  and  Kawamura,  these  experiments  reaffirmed 
their  findings  with  respect  to  rate  of  sand  transport. 

2.  The  average  flying  distance  of  the  sand  particles  was  found  to  be 
much  greater  than  that  found  by  previous  investigators,  but  the 
difference  could  pcssibly  be  due  to  the  method  of  calculation. 

3.  The  experiments  seem  t<  prove  that  sand  movement  has  «  negligible 
response  time  to  a  change  in  wind  velocity. 
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ADDENDUM  I 

STUDIES  WITH  0,30  MM  DIAMETER  SAND 


The  section  in  the  previous  study  dealing  with  the  variation  in  the 
rate  of  sand  transport  and  flying  distance  with  the  wind  velocity  has  been 
repeated  with  a  sand  of  smaller  grain  size. 

The  sand  used  in  this  new  study  had  a  mean  diameter  of  0.30  m*  and  a 
rather  narrow  range  of  grain  sizes  (Figure  1). 

The  experimental  procedure  waa  exactly  the  same  and  U*  was  again 
calculated  by  the  formula: 

U  m  6.13  U*  log  |t  +  U* 

A  sand  feed-in  was  used  in  all  runs  in  order  to  create  an  artificial 
impact  on  the  sand  bed.  This  compensates  for  the  size  of  the  bed  so  that 
the  results  obtained  approach  those  which  would  have  been  given  by  a  sand 
bed  of  infinite  length. 

The  threshold  shear  velocity  (with  sand  feed-in>  was  16  em/sec.  This 
value  agrees  fairly  well  with  the  value  obtained  from  the  Bagnold  formula 
concerning  the  impact  threshold:  _ 

U^.  *  0,08  \J  (2—— E.)  gd 
1  P 

which  gives 

IK  •  18  cm/sec  for  H  •  0.30  mm. 

Wt 

The  coordinates  of  the  focal  point  were  fou^d  to  be: 

Z'  ■  0.010  ft 

U'  -  9  ft/sec 

As  for  the  sand  used  in  original  experiment  (d  «  0.44  mm.),  thj  focal 
point  agrees  quite  well  with  Zingg's  estimate  of, 

Z*  ■  0.0099  ft 

U*  »  8.3  ft/sec 

The  rate  of  sand  transport  is  plotted  against  the  sheer  velocity  in 
Figure  2.  This  curve  intersects  the  curve  obtained  with  the  0.44  mm  sand. 


X-l 


This  fact  could  have  been  theoretically  foreseen,  since  according  to 
Bagnold  formulae; 


the  threshold  shear  velocity  ,  and  the  rate  of  sand  transport,  q,  both 
increase  with  the  grain  size  d. 

Figure  3  was  prepared  to  present  a  comparison  of  the  above  data 
(Figure  2)  with  the  experimental  results  of  Bagnold  and  Kawamura.  The 
three  curves  differ  considerably  although  the  sand  in  each  case  had 
almost  the  same  mean  diameter.  These  variations  might  be  explained  by 
differences  in  the  experimental  equipment  (length  of  the  tunnel,  presence 
or  absence  of  sand  feed-in,  and  above  all  the  length  of  the  sand  trap). 

Also  the  wind  velocity  distribution  is  sometimes  difficult  to  define 
exactly,  and  the  slightest  deviation  in  its  slope  greatly  influences  the 
calculated  value  of  U*. 

In  the  study  with  the  sand  of  0.44  mm,  diameter  the  Bagnold  and 
Kawamura  formulae  could  be  used  to  describe  the  experimental  data. 

However,  in  the  present  ease  it  was  impossible  to  find  a  constant  in 
these  formulae  which  would  permit  an  adequate  description  of  the  experi¬ 
mental  data.  Figures  4  and  5  show  the  best  description  which  could  be 
obtained.  The  explanation  of  this  fact  was  found  by  plotting  the  ex¬ 
perimental  data  on  log-log  paper  (Figure  6).  The  data  follow  along  a 
straight  line  the  slope  of  which  is  about  2.8  instead  of  3  as  theoretically 
found  by  Bagnold  and  Kawamura.  A  similar  disagreement  had  been  already 
noticed  by  Horikawa,  but  the  difference  is  too  small  in  this  case  to 
question  the  theoretical  formulae.  In  fact  the  slight  disagreement  which 
appears  mostly  for  higher  values  of  the  wind  velocity  could  be  due  to  the 
v.n». ensured  amount  of  sand  transported  in  suspension  by  the  wind:  the 

of  sand  in  suspension  probably  is  not  negligible  for  such  a  fine 
sand  at  the  high  velocities. 

Figure  7  shows  a  typical  horizontal  distribution  of  sand  in  the 
horizontal  trap.  This  curve  can  be  used  to  calculate  the  average  flying 
distance  of  the  sand  particles.  For  UA  •  40  cm/scc,  the  average  flying 
distance,  L,  equals  9  inches  which  is  much  higher  than  the  values  obtained 
by  previous  investigators  which  are  generally  about  2  inches.  No  clear 
explanation  could  be  suggested  for  this  discrepancy,  but  the  shortness  of 
the  traps  used  by  Bagnold,  Kawamura  and  Horikawa  might  be  a  source  of  error 
in  their  findings.  Also  the  possible  presence  of  large-scale  turbulence 
above  the  traps  might  invalidate  the  present  results. 
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FIGURE  I  -  GRAIN  SIZE  DISTRIBUTION 
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FIGURE  3-  RATE  OF  SAND  TRANSPORT 
Bagnold  formula: 


Bognold 
formula 
with  C  =  2.8 


U*in  cm  /sec 

FIGURE  4  -  COMPARISON 
AND 

BETWEEN  EXPERIMENTAL  RESULTS 
BAGNOLD  FORMULA 

if 


Kowamura  formula, 
q  =  k  |  (  U*-U*,)(U*  +  U.t)# 
U*t=  16  cm  /sec 


Kowomura  formula 
with  K«  3.1 


U#  in  cm/^ec 

FIGURE  5-COMPARISON  BETWEEN  EXPERIMENTAL  RESULTS  AND 

KAWAMURA  FORMULA 
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Conclusions 

This  additions!  stud"  for  the  smaller  ,;'jin  size  has  simply  confirmed 
the  conclusions  of  the  original  experiments  on  the  larger  grain  size. 

The  theoretical  aspect  of  the  subject  does  not  seem  to  pose  any 
problem,  hut  further  study  is  needed  to  explain  the  discrepancies  in  the 
experimental  data  obtained  by  the  various  investigators  concerning  the 
rate  of  sand  transport  and  the  flying  distance  of  the  particles. 


ADDENDUM  II 


SAND  TRANSPUT  RV  WIND 
STUDIES  WITH  SAND  C  <0.1  '15  MM  Hi  A  METER  ) 

by 

Abdel -Lat if  Kadib 

This  study  on  the  variation  in  the  rate  of  sand  transport  with  the 
wind  velocity  was  made  with  a  sand  having  a  mean  diameter  of  0.145  mm  and 
grain  size  distribution  as  shown  in  figure  1. 

The  experimental  procedure  essentially  was  the  same  as  that  used  by 
Delly^1)  in  tests  on  Sand  A  <D  =  0.44  mm),  and  Sand  B  (D  =  0.30  nan).  The 
only  difference  in  the  procedure  was  that  to  obtain  wind  velocities  as  low 
as  10  ft/sec,  it  was  essential  to  open  some  of  the  top  covers  of  the  wind 
tunnel.  Vertical  wind  profiles  were  measured  for  different  combinations 
of  cover  openings.  The  results  of  these  tests  are  3hown  in  Figure  2. 

Sand  was  fed  into  the  upper  end  of  the  tunnel  in  all  runs. 

The  shear  velocity  U^  was  calculated  using  Zingg’s  formula: 

U  n  6.13  U#  Log  *  U' 

The  experimental  threshold  velocity  (with  sand  feed)  was  22.0  cm/sec 
and  that  calculated  by  the  Bagnold  formula  was  17.6  cm/sec.  The  difference 
between  the  observed  and  the  calculated  values  of  the  threshold  velocity 
is  due  to  the  use  of  a  constant  value  for  A  in  the  Bagnold  formula: 

u*  =  A  Z7r~S'  6d 

t  p 

Bagnold  showed  that  the  value  of  A,  for  a  grain  size  of  0,20  mm  or  less, 
is  not  a  constant.  The  coordinates  of  the  focal  point  of  the  vertical 
velocity  distribution  curves  were  obtained  from  Figure  3. 

Z'  »  0.0125  ft 
U'  *  6.4  ft/sec 

These  values  did  not  agree  with  that  obtained  from  the  Zingg  equation,  i.e. 

Z'  =  0.0048  ft 
U'  =  4.25  ft/sec 

Table  1  gives  the  experimental  data  obtained  from  these  tests.  These 
data  are  plotted  in  Figure  4. 


(1)  Belly,  Pierre-Yvcs,  Sand  Movement  by  Wind,  Univ.  of  Calif..  Institute 
of  Engrg  Res.,  Series  77,  Issue  7,  July  1962,  (included  as  main  text 
and  Addenda  I  8i  III  to  this  technical  memorandum). 
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Figure  5  shows  the  experimental  values  of  sand  transport  compared  with 
values  calculated  by  the  Bagnold  and  K  iwamura  formulas.  In  this  comparison, 
a  value  of  C  —  1.5  was  used  in  ihe  Bagnolu  formula.  This  is  the  lower  value 
for  the  value  of  C  recommended  by  Bagnold  and  applies  to  a  uniform  grain 
size.  For  the  Kawamura  formula  a  value  of  K.  =  1.0  had  to  be  selected  in 
order  that  the  calculated  curve  would  pass  near  the  experiential  points. 

Comparison  with  O'Brien  and  Rindlaub  formula 

(2) 

O'Brien  and  Rindlaub  '  proposed  the  following  formula  from  data  de¬ 
rived  from  the  field  tests: 

G  =  0.036  (for  Us  >  20  ft/sec) 

where  G  is  the  rate  of  movement  of  dry  sand  in  pounds  per  day  passing  an 
imaginary  line  1  foot  in  length  drawn  perpendicular  to  the  wind,  and  II 
is  the  wind  velocity  5  ft,  above  the  sand  surface  in  ft/sec.  During  tne 
field  tests  the  sand  caught  in  the  sand  trap  had  a  median  diameter  between 
0.17  mm  and  22  mm. 

In  order  to  compare  the  test  results  with  the  O'Brien  and  Rindlaub 
formula,  the  velocities  5  ft.  above  the  sand  bed  in  the  present  tests  were 
calculated  using  Zingg's  formula,  i.e., 

U  =  6.13  U*  Log  |,  +  U' 

This  calculation  was  made  for  Sand  A  (D  *  0.44  mm).  Sand  B  (D  =  0.30  mm) 
and  Sand  C  (D  =  0.145  mm).  Tables  2,  3  and  4  show  the  calculations  for 
Sands  A,  B,  and  C,  respectively.  Figure  6  shows  the  plots  of  the  experi¬ 
mental  data  compared  with  the  O'Brien  and  Rindlaub  formula. 

The  experimental  values  for  the  sand  transport  are  compared  with  the 
O'Brien  and  Rindlaub  formula  as  follows: 


Sand  A  (D  =  0.44  mm).  The  data  follow  a  straight  line  the  slope  of  which 
is  about  6.00  instead  of  3.00  as  found  by  O'Brien  and  Rindlaub  (Figure  6a). 
Also  the  value  of  the  constant  (0.036)  in  their  formula  was  found  to  be 
about  0.76  x  1CT6  for  sand  A,  This  limits  the  use  of  their  equation  to 
sand  with  a  grain  diameter  similar  to  what  they  tested  (D  =  0.165  nun). 

For  Sand  A  (D  =  0.44  mm)  the  best  description  is 


G  =  0.76  x  10”6  U5b‘°° 


Sand  B  (P  =  0.3  mm).  Experimental  data  follow  a  straight  line  parallel  to 
the  O'Britn  and  Rindlaub  curve  but  shifted  to  the  left  (Figure  6b),  The 


(2)  O'Brien,  M.  P.  and  Rindlaub,  B.D. ,  The  transportation  of  sand  by 
wind,  Civil  Engrg.,  May  1936,  pp  325-326. 
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value  of  the  constant  for  Sand  B  was  found  to  he  0.065.  For  Sand  B 
(D  =  0.3  mm) ,  the  best  description  is 

G  =  0.065  x  U53‘°° 


Sand  C  (D  =  0.145  mm).  Figure  6c  shows  a  scatter  of  the  experimental  data 
around  the  O'Brien  and  Rindlaub  curve.  The  scatter  does  not  seem  to  be 
significant.  This  Sand  was  the  only  one  which  agrees  with  the  0”Brien  and 
Rindlaub  formula.  This  agreement  probably  is  because  Sand  C  (  D  =  0.145  mm) 
has  approximately  the  same  order  of  magnitude  of  grain  diameter  as  found  in 
the  O'Brien  and  Rindiaub  tests  (D  =  0.105  mm). 


CONCLUSION 


1.  The  value  of  threshold  velocity  is  best  determined  by  experiment 
when  the  grain  size  is  less  than  0.20  mm. 

2.  The  value  of  C  in  Bagnold’s  formula  has  less  scatter  than  the  co¬ 
efficient  K.  in  Kawamura's  formula. 

3.  The  use  of  the  O'Brien  and  Rindlaub  formula  must  be  limited  to 
sand  having  the  same  grain  diameter  of  that  tested. 


Table  l 


;  velocity 

Shear  velocity 

Rate  of 

«=  1.0  ft. 

in  cm/sec 

transport 

i  ft/sec 

in  gr/cm-sec 

14.70 

24.40 

0.00167 

16,00 

27.80 

0.0277 

18.20 

33.60 

0.0590 

19.40 

36.80 

0.0880 

21.00 

41.00 

0.114 

26.60 

55.50 

0.278 

28.00 

59.50 

0.189 

29.20 

62.60 

0.270 

30.00 

64.50 

0.365 

33.20 

73.00 

0.475 

34.20 

76.00 

0.530 

35,40 

79.00 

0,645 

36.00 

80.50 

0.736 

37.50 

84.50 

0.950 

38.90 

88.00 

0.945 

H-3 


Table  2 


Sand  A 

(D  =  0,44  mm) 

U'  = 

13  ft/sec 

Z'  = 

0.0135 

U5  = 

U*  6'13  Lo«  0.0135  +  13 

,0 

u5  « 

15.7  U#  +  13.00 

«* 

15.7  U* 

«5 

ft/sec 

q 

1 b/f t-day 

1.48 

23.2 

36.2 

1820.0 

1.48 

23.2 

36.2 

1690.0 

1.52 

23.8 

36.8 

2220.0 

1.65 

25.8 

38.8 

2920.0 

1.68 

26.4 

39.4 

3160.0 

1.71 

27.0 

40.0 

3360.0 

1.78 

28.0 

41.0 

4050.0 

1.81 

28.4 

41.4 

4540.0 

1.91 

30.0 

43.0 

5260.0 

2.1 

33.0 

46.0 

5850.0 

0.99 

15.5 

28.5 

69.5 

1.15 

18.1 

31.1 

610.0 

1.25 

19.6 

32.6 

1055.0 

1.28 

20.2 

33.2 

1280.0 

1.35 

21.2 

34.2 

1350.0 

1.52 

23.8 

36.8 

1625.0 

1.65 

25.8 

38.8 

2930.0 

Table  3 

Calculations  for  U5  ft  using  Zingg'a  Formula 
Sand  BCD-  0.30  mm) 

U'  =9  ft/sec 
Z'  =  0.01C 

U5  =  UA  6.13  Log10  ^  ♦  9.00 

=  16.5  U*  +  9 


u* 

ft/sec 

16.5  U+ 
ft/sec 

"5 

ft/sec 

q 

lb/ft -day 

0.545 

8.70 

17.70 

0 

0.825 

13.60 

22.60 

665.0 

0.99 

16.30 

25.30 

930.0 

1.085 

17.80 

26.80 

1270.0 

1.180 

19.50 

28.50 

1680.0 

1.315 

21,70 

30.70 

1680.0 

1.35 

22.30 

31.30 

2140.0 

1.42 

23,40 

32.40 

2320,0 

1.58 

26.00 

35.m 

2960.0 

1.90 

31,40 

40.40 

4400.0 

2.10 

34.80 

43,80 

5600,0 

Table  4 


Calculation  of  ^  using  the  Logarithmic  Formula 
Sand  C  (D  =  0.145  mm) 

Using  the  Modified  Wind  Velocity  Distribution 
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FIGURE  2.  EFFECT  OF  TOP  OPENINGS  ON  VERTICAL  WIND  PROFILE 
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Wind  Velocity  5  Feet  Above  Send  Bed,  ft /sec 

FIGURE  6.  COMPARISON  BETWEEN  EXPERIMENTAL  DATA  AND 
THE  O’BRIEN  ANO  RINDLAUL8  FORMULA 


ADDENDUM  1X1 

INFLUENCE  OP  MOISTURE  OK  THE  THRESHOLD  OF  SAND  MOVEMENT 


INTRODUCTION 

When  a  flow  of  air  over  a  flat  bed  of  loose  grains  is  gradually  in¬ 
creased,  there  comes  a  certain  moment  when  the  sand  grains  are  (Hit  into 
motion  by  the  force  of  the  wind  (which  can  be  measured  in  terms  of  shear 
velocity  U*).  This  critical  value  of  the  wind  velocity  is  cslled  the 
threshold  velocity,  snd  the  corresponding  vslue  of  the  shear  velocity,  the 
threshold  shear  velocity  is  U** .  This  threshold  velocity  depends  mainly 
upon  the  characteristics  of  the  sand  and  of  its  surface. 

This  subject  has  been  investigated  by  many  authors,  among  them, 
Jorissen^^Jeffreys'2}  Chepil'3)(  zingg^*',  but  chiefly  by  Bagnold”’8'9^ 
in  the  particular  case  of  sand. 

Fluid  and  Impact  Threshold 

The  fluid  threshold  is  the  critical  value  of  the  wind  velocity  which 
has  to  be  reached  in  order  to  initiate  movement  in  the  hypothetical  case 
of  an  extremely  flat  sand  surface.  The  movement  of  the  sand  grains  in 
this  case  is  caused  only  by  the  drag  of  the  wind. 

The  impact  threshold,  on  the  other  hand,  occurs  when  the  sand  is  kept 
disturbed  by  the  "impact"  of  oncoming  grains  upon  it.  It  is  the  impact 
threshold  which  is  generally  observed  under  natural  conditions,  for  there 
is  always  a  temporary  stronger  wind  or  some  irregularities  of  the  bed 
(pebbles,  small  mounds)  which  disturb  the  flow  and  increase  the  velocity 
locally.  Both  cause  a  local  or  temporary  sand  transport.  The  impact  of 
grains  already  in  motion  on  immobile  sand  helps  the  wind  drag  to  put  the 
latter  into  motion;  and  this  process  repeats  itself.  After  a  short  time, 
the  sand  movement  extends  over  the  entire  surface.  Thus,  a  continuous 
saltation  of  grains  can  be  maintained  for  an  indefinite  distance  downstream 
by  a  wind  of  feebler  strength  than  the  fluid  threshold.  The  threshold  wind 
strength  at  which  this  occurs  marks  the  critical  stage  at  which  the  energy 
supplied  to  the  saltating  grains  by  the  wind  just  balances  the  energy 
losses  due  to  friction  when  grains  strike  the  bed. 

Impact  and  Fluid  Threshold  for  Uniform  Sand 

As  a  large-scale  phenomenon,  the  wind  over  an  open  dune,  and  the  air 
flow  in  a  wind  tunnel  (at  a  distance  from  the  tunnel  entrance  for  ar>  open- 
circuit  type  tunnel)  can  be  considered  an  fully  turbulent.  Dealing  now 


*For  references  cited  in  Addendum  III  see  page  III-24 
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with  small-scale  flow  over  and  around  the  individual  grains  on  the  surface, 
the  Reynolds  number  characterizing  the  flow  is 

U*d 

v 

where  U*  is  the  shear  velocity,  d  the  mean  size  of  the  sutface  roughness 
(which  is  of  the  same  order  as  the  grain  diameter),  and  v  the  kinematic 
viscosity  of  the  fluid  <0,14  for  air  in  c.g.s.  units). 

If  Hi-5  >  3.5,  the  surface  is  rough  and  the  threshold  shear  velocity 

which  must Vbe  attained  by  the  fluid  before  it  can  move  any  surface  grains 
varies,  as  shown  by  Bagnold,  as  the  square  root  of  the  grain  diameter, 
according  to  the  equation, 

U*t  ’  At/iVe  fid 

where  a  is  the  density  of  the  grain  material  and  p  is  the  density  of  air. 

A  is  a  coefficient  which  equals  0.1  in  the  case  of  the  fluid  threshold. 

Since  the  impact  threshold  is  lower  than  the  fluid  threshold,  the  coeffi¬ 
cient  A  will  be  slightly  smaller,  and  approximates  the  value  of  0.08. 

For  smaller  grains,  when  V*~<  3.5,  the  surface  approaches  the 
smooth  condition  and  the  coefficient  A  ia  no  longer  constant  but  increases 
as  the  grains  become  smaller  and  smaller.  Figure  1  gives  the  variations 
of  Fluid  and  Impact  Threshold  with  grain  size  as  found  by  Bagnold. 

The  Fluid  Threshold  for  Natural  Sand 

Usually  the  natural  sand  is  a  mixture  in  which  one  size  of  grain  pre¬ 
dominates,  and  in  which  the  proportions  by  weight  of  grains  of  greater  and 
smaller  diameter  decrease  ss  the  size  departs  from  that  of  the  predominant 
grains.  If  the  material  ia  well  mixed  and  is  spread  out  over  the  ground, 
the  surface  may  be  assumed  to  contain  exposed  grains  in  the  proportions  in 
which  they  exist  in  the  body  of  the  material.  Most  of  the  fine  grains  lie 
in  crevices  between  the  larger  ones,  and  are  screened  by  these  from  the 
drag  of  the  wind.  Apart  from  a  very  temporary  movement  on  the  part  of  the 
few  most  exposed  fine  grains  for  which  the  threshold  velocity  is  smsllcr 
than  for  the  mean  diameter,  the  initial  threshold  wind  velocity  is  that 
corresponding  to  the  predominant  diameter. 

If  the  wind  is  not  increased  above  the  initial  fluid  threshold  strength, 
sand  movement  goes  on  until  the  grains  of  the  perdominant  and  smaller  sizes 
have  all  been  carried  away  from  the  exposed  surface,  leaving  only  those  of 
larger  diameter.  Then  the  motion  ceases.  By  raising  the  wind  velocity,  a 
further  temporary  movement  is  produced,  and  so  on.  Finally  the  result  is 
a  sand  bed  covered  by  a  surface  layer  containing  all  the  largest  grains 
which  were  present  in  the  removed  layers.  If  vTie  wind  strength  is  raised 
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FIGURE  I-  VARIATION  OF  THE  THRESHOLD 
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still  further  until  the  largest  grains  begin  to  move,  the  motion  is  no 
longer  temporary,  but  goes  on  indefinitely.  This  ultimate  threshold  is 
tnar  corresponding  to  ihe  largest  grains  lut-aeiil  in  the  bed. 

Since  the  wind  may  drop  before  the  ultimate  threshold  is  reached, 
the  sand  may  be  left  in  any  state  of  surface  arrangement  and  the  wind 
strength  required  to  move  it  again  may  be  anything  from  the  initial  to 
the  ultimate  threshold,  tor  normal  sand,  however,  the  biggest  grains 
oresent  are  usually  not  more  than  twice  the  predominant  diameter,  so 
that  the  ultimate  threshold  is  commonly  exceeded,  and  the  transient 
stages  occur  only  over  a  small  range  of  wind  velocity. 

The  Impact  Threshold  for  Natural  Sand 

Physically  the  impact  threshold  marks  the  critical  stage  at  which 
the  energy  supplied  to  the  saltating  grains  begins  to  balance  the  energy 
losses  due  to  friction  when  the  grains  hit  the  bed.  There  appears  there¬ 
fore  to  be  little  direct  connection  between  the  impact  threshold  and  the 
conditions  at  the  surface  which  determine  the  fluid  threshold.  Bagnold 
found  that  for  sand  of  mixed  grain  size  the  impact  threshold  is  approxi¬ 
mately  that  corresponding  to  the  predominant  diameter  of  the  grains  in 
surface  creep. 


PURPOSE  OF  THE  EXPERIMENT 

Measurements  on  sand  movement,  and  particularly  of  the  threshold 
velocity  show  a  great  scatter  for  both  field  and  laboratory  experiments. 
Wind  velocity  and  grain  size  of  the  sand  are  certainly  the  major  factors 
of  sand  movement ,  but  due  to  the  inherent  complexity  of  the  problem 
many  other  variables  can  play  minor  roles  and  could  be  responsible  for 
the  observed  scatter.  Below  are  summarized  the  possible  factors; 

-  related  to  the  wind  Temperature 

Humidity 

-  related  to  the  sand  Structure 

Texture 

Moisture -Content 

-  related  to  the  surface  Roughness 

Settlement 

Temperature 

The  moisture  adsorbed  by  the  sand  holds  its  grains  together. 
Evidence  for  this  esn  be  obtained  from  the  different  angles  of  repose 
of  dry  and  wet  sand.  It  is  therefore  easy  to  foresee  that  the  moisture 
renders  the  extraction  of  the  grains  of  a  sand  bed  more  difficult.  The 
purpese  of  this  experiment  is  to  evaluate  the  general  magnitude  of  this 
effect. 
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Ip,  order  to  study  the  thteshold  of  sand  movement  for  different  sir 
humidities  it  is  preferable  to  use  a  closed-circuit  type  wind-tunnel. 
Contrary  to  the  open-circuit  type,  such  •  tunnel  does  not,  at  least  in 
principle,  require  a  constant  addition  of  water  vapor. 

In  order  to  simulate  full-scale  sand  movement  in  a  tunnel  it  is 
necessary  that  the  air  stream  have  a  velocity  magnitude  and  distribution 
and  turbulence  structure  similar  to  those  of  natural  wind  near  the 
ground,  and  that  the  working  section  be  large  enough  to  minimize  wall- 
interference. 

For  the  particular  study  of  the  fluid  threshold  of  sand  movement 
the  bed  should  be  limited  in  length,  and  no  sand  enter  from  above  the 
up-wind  extremity  of  the  bed,  because  only  in  such  a  condition  can  the 
true  fluid  threshold  be  examined  as  distinct  from  the  impact  threshold. 

A  wind  tunnel  located  in  Building  160  at  the  University  of  Cali¬ 
fornia  Richmond  Field  Station  fulfilled  the  above  conditions  and  was 
used  in  the  present  experiment.  The  elevation  and  plan  of  this  tunnel 
is  shown  in  Figure  2.  The  working  section  is  on  the  suction  side  of 
the  fan  and  is  shown  in  Figure  3  with  a  sand  bed  ready  for  testing. 

This  chamber  is  1  foot  high,  2  feet  wide  and  3  feet  long.  One  side  was 
n«;de  of  plastic  to  permit  visual  observations  during  testa. 

The  motor  operating  the  fan  was  always  set  at  a  constant  rotational 
epeed.  A  wind  speed  range  of  17  to  43  ft/sec  could  be  obtained  by  con¬ 
trolling  a  sliding  valve  placed  on  the  exhaust  side  of  the  fan. 

The  sand  was  placed  in  a  box,  2  inches  deep  and  23  inches  long, 
which  could  be  removed  from  the  working  chamber.  The  surface  of  the 
sand  was  brought  to  a  height  equal  to  that  of  the  upwind  and  downwind 
floor  of  the  tunnel. 

Since  the  experiments  were  never  conducted  for  high  rates  of  sand 
transport,  most  of  the  blown-off  sand  could  be  collected  in  a  sand  trap, 
2  inches  long,  placed  immediately  downwind  of  the  end  of  the  sand  bed. 
The  sand  which  was  not  caught  in  this  trap  was  deposited  in  the  vicinity 
of  the  valve,  so  tuat  tiii-  air  was  free  of  sand  particles  when  arriving 
again  in  the  working  section. 


Wind  Velocities  Measurements 

It  is  essential  when  measuring  the  velocity  gradient  at  which  the 
fluid  begins  to  cause  sand  movement  that  a  sufficient  length  of  sand 
surface  exist  upwind  of  the  point  of  measurement,  so  that  ;he  fluid  may 


FIGURE  2  -  WIND  TUNNEL 


FIGURE  3-  THE  WORKING  CHAMBER.  THE  SAND  BOX 

AND  ITS  TRAP  ARE  SHOWN  BELOW,  THE  VELOCITY 
MEASURING  EQUIPMENT,  WET  AND  DRY  THERMO¬ 
METER  AND  THE  RHEOSTAT  WHICH  CONTROLS  THE 
HEATERS  ARE  ABOVE. 
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have  an  opportunity  to  adjust  the  velocity  distribution  appropriate  to 
the  texture  of  the  sand  bed.  Otherwise,  except  at  levels  immeasurably 
Close  TO  the  surface,  the  velocity  uisl  i  ibui.iajw  will  l»c  ilia l 
smooth  floor  further  upwind. 

Wind  velocities  were  measured  with  a  standard  Prandtl  type  pitot 
tube  which  was  attached  to  a  point  sage  and  introduced  into  the  air 
stream  through  the  top  of  the  flume.  In  order  that  the  velocity  dis¬ 
tribution  adjust  to  the  sand  bed,  the  pitot  tube  was  placed  along  a 
vertical  which  was  almost  at  the  end  of  the  sand  bed.  The  pitot  tube 
was  connected  to  an  Ellison  type  draft  gage  having  a  range  of  one  inch 
of  water  and  graduated  into  divisions  of  0.01  inch.  The  pitot  tube  i’.ad 
a  "coefficient"  of  1.  The  wind  velocities  therefore  were  calculated 
by  the  formula,  _ 

»•  J 

where  p  is  the  differential  pressure  at  the  pitot  tube,  and  p  the 
density  of  air,  and  g  the  acceleration  of  gravity.  Atmospheric  pressure, 
air  humidity  and  air  temperature  were  taker  into  account  in  determining 
the  air  density. 

Sand  Dampening 

In  order'  to  study  the  movement  of  moist  sand  by  wind,  during  a 
preliminary  trial  the  sand  was  dampened  by  spraying  some  wmter  on  its 
surface  and  a  wind  of  normal  humidity  was  allowed  tc  blow  through  the 
tunnel.  It  was  immediately  noted  that  this  method  was  inappropriate 
because  tne  wind  dried  the  surface  layer  of  the  sand  bed  invalidating 
the  results.  Thus  the  following  procedure  was  adopted. 

1)  Study  the  movement  of  sand  by  winds  of  various  humidities. 
The  wind  was  allowed  to  blow  over  the  sand  for  a  sufficiently  long  time 
to  dampen  the  sand  surface, 

2)  Study  the  movement  of  sand  by  wind  for  higher  water  content 
of  the  sand,  by  adding  water  directly  to  the  sand.  In  this  case  it  was 
necessary  to  have  the  wind  saturated  with  water,  so  that  the  surface 
layer  of  the  sand  does  not  dry  quickly. 

There  are  several  methods  of  measuring  the  moisture  content  cf  a 
soil.  Among  them  are  electrical  methods  (generally  based  on  resistivity 
measurements)  and  radioactive  methods.  They  all  require  more  or  less 
complicated  equipment ,  and  they  all  have  to  be  calibrated  with  reference 
to  the  direct  method  consisting  of  weighing  samples  before  and  after 
drying.  They  are  therefore  less  accurate  than  the  direct  method.  The 
direct  method  is  generally  not  feasible  b-cause  of  the  time  element 
involved.  As  an  oven  was  made  available  by  the  Sanitary  Engineering 
Laboratory  of  the  University,  and  since  time  was  not  a  controlling 
factor,  the  direct  method  was  therefore  used. 
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Sampling  was  done  inside  the  tunnel  and  sample  boxes  were  immedi¬ 
ately  sealed  (also  inside  the  tunnel  in  order  to  retain  in  the  box  the 
air  of  the  tunnel).  The  samples  were  then  taxes  to  tne  Sanitary  Engin¬ 
eering  Laboratory,  weighed  (using  a  balance  sensitive  tc  a  1/100  of  a 
gram)  and  then  placed  in  an  oven  at  105°  Centigrade.  Twenty-four  hours 

later  the  sample  containers  were  again  weighed  and  the  water  content,  w, 
calculated  as, 

■  in  %  .  Wet  weight  -  dry  weight  1QQ 
dry  weight 


Humidification  of  Air 

Water  vapor  was  added  to  the  air  inside  the  tunnel  by  using  a  pan 
of  water  heated  by  three  submerged  666-wmtt  electric  heaters.  The  free 
surface  of  the  water  had  an  area  of  about  3  ft2.  While  heating  the 
water,  the  fan  was  used  from  time  to  time  to  help  circulate  the  air  in 
order  to  obtain  a  homogeneous  mixture  of  water  vapor  and  air. 

The  air  humidity  is  defined  as, 

m  amount  of  water  vapor  present 

amount  of  water  vapor  for  saturation 

and  was  determined  by  using  a  wet  and  dry  bulb  thermometer  and  tables 
published  by  the  U.  S.  Weather  Bureau(5).  The  apparatus  was  olaced  in 
the  working  section  and  the  wet  bulb  kept  constantly  wet  by  a  wick 
dipped  in  water. 

It  was  very  difficult  to  keen  the  air  at  a  particular  humidity 
during  tests,  because  the  water  vaoor  condensed  on  the  wall  of  the 
tunnel  and  around  the  fan  when  air  is  circulating  in  the  tunnel.  Thus, 
as  much  as  oossible,  the  experiments  were  conducted  with  air  at  its 
normal  humidity.  The  natural  variations  in  air  humidity  (45  to  70%) 
were  sufficiently  large  to  give  significant  results.  The  artificial 
humidification  was  used  to  experiment  with  air  of  humidity  approaching 
saturation  and  in  this  case  the  air  humidity  was  taken  as  the  average 
between  the  measurements  made  at  the  beginning  and  at  the  end  of  each 
run. 

Measurement  of  the  Threshold 

A  natural  sand  having  a  mean  diameter  of  0.44  mm  was  usee  ’  this 
study.  This  sand  had  a  relatively  wide  range  of  grain  sizes  as  can  be 
seen  in  Figure  4.  As  already  noted  the  threshold  will  therefore  occupy 
^  relatively  wide  range  of  wind  velocities.  It  was  decided  to  take 
account  only  of  the  "ultimate  threshold",  i.e.,  the  lowest  wind  strength 
which  gives  a  general  movement  of  the  sand.  This  choice  was  dictated 
by  the  following  considerations: 
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FIGURE  4 -GRAIN  SIZE  DISTRIBUTION 


(a)  the  ultimate  threshold  has  more  significance  than  the 
initial  threshold  as  far  as  natural  conditions  are  concerned.  It  is 
the  lowest  value  of  the  wind  strength  which  given  a  measurable  sand- 
removal  . 


(b)  the  value  of  the  ultimate  threshold  does  not  deoend  on  the 
past  history  of  the  sand  bed  as  is  the  case  for  the  initial  threshold, 
since  for  the  former,  grains  of  all  diameters  can  be  moved  by  the  wind. 

(c)  the  wind  velocity  measurements  are  easier  to  make  because 
the  state  of  ultimate  threshold  lasts  indefinitely,  whereas  ti.c  initial 
threshold  lasts  only  as  long  as  the  grains  of  the  predominant  diameter 
are  present  on  the  surface  layer. 

The  accurate  definition  of  the  ultimate  threshold  is,  however,  a 
difficult  matter.  Thus,  instead  of  a  visual  observation,  the  threshold 
was  defined  as  occurring  when  the  rate  of  sand  caught  in  a  small  trap 
beyond  the  end  of  the  sand  bed  had  a  certain  value  (0.04  gram/c»-sec). 
When  this  rate  had  not  been  obtained  directly,  an  interpolation  or 
extrapolation  method  was  used. 


EXPERIMENTAL  RESULTS  AND  DISCUSS ION 


Preliminary  Remarks 

Besides  supplying  water  to  the  sand,  air  humidity  plays  a  role  in 
the  threshold  problem  by  changing  the  ratio  of  the  densities  of  sand 
particles  and  fluid.  The  Bagnoid  formula  provides  for  this  change  in 
air  density,  that  is, 

u*t  -  A  7  ”7”  *  d  • 

For  a  temperature  of  70°F  and  a  pressure  of  30  inches  of  mercury,  the 
completely  dry  air  (h  =  05k)  has  a  density  of  0.0755  lb/ft3.  For  the 
same  temperature  and  pressure,  the  air  saturated  with  water  (h  *  100%) 
has  a  density  of  0.0743.  a  being  much  bigger  than  p,  the  above  formula 
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can  be  rewritten  as. 


thus. 
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and  for  the  above  variation  of  n .  the  relative  variation  of  U, ,  is, 

H  •  *t 


d  It. 


*t 
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This  variation  is  so  small  in  comparison  with  that  which  will  be  found 
for  actual  variations,  that  the  influence  of  humidity  on  air  density  c'n 
be  neglected  with  respect  to  the  role  of  sand  moisture.  Moist  sand  found 
along  coasts  can  be  classified  into  two  categories  depending  upon  the 
origin  of  the  moisture: 


1)  Moist  sand  which  has  collected  moisture  from  the  atmospnere. 
Unlike  very  fine  particles  such  as  dust  or  loess  the  sand  does  not  rtaiily 
absorb  moisture.  As  proved  later,  there  exists  a  correlation  between  the 
water  content  of  the  sand  and  air  humidity.  Allowing  for  this  fact,  and 
since  in  the  field  it  i*  easier  to  measure  the  air  humidity  than  the  sand 
water-content,  the  air  humidity  instead  of  the  water-content  of  the  sand 
has  been  taken  as  the  variable  in  this  study. 


2 )  Moist  sand  whose  water  comes  from  sources  other  than  air 
humidity,  such  as:  rain,  rising  of  underground  water,  and  sea  water  re¬ 
maining  in  the  sand  by  wave  or  tide  action.  When  a  wind  not  saturated 
with  vapor  blows  over  such  a  sand,  it  gradually  dries  out  the  surface 
layer  of  the  sand  bed,  until  an  equilibrium  is  reached  between  humidity 
of  the  wind  and  the  water  content  of  the  surface  layer  of  the  ground. 

The  study  of  sand  movement  ir.  this  case,  where  water  contents  of  air  and 
sand,  as  well  as  wind  duration,  come  into  play,  is  very  complex.  The 
first  step  in  the  solution  of  this  problem  consisted  of  studying  the 
particular  case  of  saturated  wind  blowing  over  a  bed  of  moist  sand  (the 
water  content  of  the  sand  being  greater  than  that  which  could  be  obtained 
from  air  humidity  alone)  for  only  in  this  instance  did  the  wind  have  no 
tendency  to  dry  out  the  sand. 

Therefore,  the  initial  studies  were  on  the  threshold  of  sand  movement 
for  various  air  humidities  and  subsequently  the  threshold  of  sand  movement 
when  directly  dampened  with  water,  by  a  saturated  wind.  Finally  an  attempt 
at  generalization  was  made  by  dropping  the  air  humidity  variable,  and  a 
relationship  between  the  threshold  and  water-content  of  sand  was  developed. 

Variation  of  the  Threshold  with  Air  Humidity 


From  a  series  of  velocity  distribution  curves  as  given  in  Figure  5, 
was  determined  for  different  mean  velocities  and  different  air  humid¬ 
ities.  Figure  6  shows  the  relationship  between  the  shear  velocity,  Ut, 
and  the  mean  velocity,  U.  This  relationship  is  quite  consistent  and  can 
be  expressed  by  a  straight  line.  Thus,  it  can  be  said  that  the  wind  drag 
is  practically  uninfluenced  by  air  hlm^^dil.;•.  In  th*  later  tests,  only 
the  value  of  the  mean  velocity  was  recorded  and  the  value  of  was 
calculated  from  this  graph. 


m-!3 


IE  -  14 


20 


FIGURE  6 


25  30 

Mean  Velocity  U  in  ft/sec 


-  RELATIONSHIP  BETWEEN  U  AND  U*  FOR  DIFFERENT 
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Note*.  The  shear  velocity  U,  was  calculated  from  the  Prondtl 
equation  with  k  *  0.401  Ux  *  5.75  U*  loq  ^ 
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Since  the  wind  velocity  never  exceeded  the  threshold  value  by  an 
appreciable  amount,  it  was  very  difficult  to  determine  the  elevation  of 
the  focal  point  of  the  velocity  distributions  (they  are  almost  parallel), 
hut  its  abscissa  is  probably  about  20  ft/sec.  This  is  scmewhat  more  than 
the  value  of  13  ft/sec  predicted  from  Zingg's  formula1^)  The  Zingg 
formula  was  well  verified  during  previous  experiments  using  the  same  sand 
but  in  a  much  longer  tunnel*,  so  that  the  present  difference  may  be  due 
to  the  shortness  of  the  bed  which  does  not  allow  the  wind  to  reach  an 
equilibrium  profile.  The  water-content  of  sand  obtained  during  the  vari¬ 
ous  runs  is  plotted  against  the  corresponding  air  humidity  in  Figure  7. 

The  scatter  is  important,  but  there  is  a  general  tendency  for  the  points 
to  follow  along  a  straight  line  (which  obviously  should  pass  through 
the  origin). 

As  explained  above,  the  threshold  velocity  was  found  by  investigat¬ 
ing  the  initial  stage  of  curves  for  the  rate  of  sand  transport.  The  three 
sample  curves  in  Figure  8  clearly  indicate  the  change  in  the  rate  of  sand 
transport  with  humidity.  From  such  curves  the  data  on  the  variation  of 
the  threshold  shear  velocity  distinctly  show  an  increase  with  air  humidity, 
and  the  relationship  is  nearly  linear  (Figure  9).  The  few  points  which 
correspond  to  a  threshold  value  manifestly  lower  than  the  average  were 
probably  obtained  during  runs  in  which  the  absorbtion  of  moisture  by  the 
sand  had  not  reached  its  equilibrium  value. 

The  increase  of  the  threshold  shear  velocity  with  air  humidity  can 
be  represented  for  the  parcicular  sand  studied,  by  the  equation 

*  28  (1  +|  h/100)  cm/sec 

where  28  is  the  value  given  by  the  Bagnold  formula  for  the  threshold 
velocity  in  cm/sec ,  for  the  sand  under  investigation.  Assuming  a  similar 
influence  of  humidity  for  sand  of  different  grain  siaes,  the  Bagnold 
formula  can  be  modified  as  follows  to  provide  for  the  air  humidity  factor. 

u*t  ■ A  u  *  I  Too 5  v-~r 

where  A  is  a  coefficient  whose  value  approximates  0.1  and  h  is  expressed 
in  percent. 

Threshold  for  Sand  of  High  Water-content 


The  maximum  amount  of  water  that  can  be  imparted  to  sand  by  atmos¬ 
pheric  humidity  is  0.25%,  For  a  water  content  greater  than  0.25%,  the 
relationship  between  U  and  U*  remains  linear  and  seems,  quite  naturally, 
to  be  a  prolongation  of  the  previous  curve  (Figure  10;. 


*  Ref. 
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FIGURE  7-  RELATIONSHIP  BETWEEN  WATER  CONTENT 
OF  THE  SAND  AND  AIR  HUMIDITY 
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FIGURE  9-  THRESHOLD  VELOCITY  AS  A  FUNCTION  OF 
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FIGURE  10-  RELATIONSHIP  BETWEEN  U  AND  U* 
(FOR  HIGH  WATER  CONTENT) 


During  these  runs  the  winrt  was  constantly  saturated  and  the  water 
content  of  the  sand,  w,  is  now  taken  as  the  variable  for  the  study  of  the 
♦hIeshold  velocity.  The  results  are  summarized  in  Figure  11.  For  a  high 
water-content  (w  >  1%}  the  wind  strength  necessary  to  initiate  sand  move¬ 
ment  becomes  more  and  more  important.  This  increase  couid  be  explained 
by  the  fact  that  the  sand  surface  becomes  very  smooth  under  wind  action, 
the  water  contained  in  the  sand  fills  up  the  interstices  between  the 
grains  making  the  extraction  of  the  grains  by  the  wind  much  mere  difficult. 
The  experiment  could  not  be  pursued  for  water-content  higher  than  4%  be¬ 
cause  the  wind  strength  necessary  to  initiate  the  movement  could  not  be 
obtained  with  the  equipment  available;  however,  one  would  ey.pt:*  that 
the  wind  strength  would  increase  very  rapidly  with  an  increase  of  water- 
content.  It  is  even  probable  that  for  the  very  high  water  content  (flooded 
sand)  the  problem  changes  aspect  and  becomes  closer  to  the  problem  of  an 
interface  between  two  fluids. 

Using  the  relationship  between  water-content  and  air  humidity  in 
Figure  7,  it  is  possible  to  complete  the  results  of  the  study  of  low  water- 
content  and  thus  find  the  relationship  between  U*t  and  w  for  the  total 
range  of  water  content  (0  to  4%).  Figure  12  shows  the  curve  obtained. 

As  one  would  have  expected,  since  the  air  humidity  in  itself  does  not  play 
a  role  in  the  sand  movement,  there  is  no  break  at  the  point  which  joins 
the  two  sets  of  data.  The  curve  suggesting  an  exponential  function  has 
been  replotted  on  semi- log  paper  in  Figure  13.  The  data  in  these  new 
coordinates  appears  to  follow  a  straight  line,  thus  indicating  that  the 
relationship  between  U*t  and  m  can  be  put  into  the  form, 

u*t  "  4  iogio  w  *  b 

where  a  and  b  are  two  constants  obtained  from  the  graph.  It  was  found 
that 

U*t  •  17  log1Q  w  +  51  cm/sec 
or , 

«  28  (0.6  log10  w  +  1-8^  cm/sec 

where  28  is  the  value  of  U.  in  cm/sec  given  by  tha  Bagnold  formula. 

These  formulae  or  the  graphs  in  Figures  J2  and  13  summarize  the 
results  obtained  in  thi3  study.  Assuming  that  moisture  affects  the 
movement  of  sand  of  different  grain  sizes  in  the  same  manner,  the  Bagnold 
formula  for  the  threshold  shear  velocity  may  be  modified  as  follows, 

U*t  -  A  j-y9-  gd  Cl. 8  ♦  0.6  log1Q  w) 
where  A  approximates  0.1  and  w  is  expressed  in  percent. 
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FIGURE  12-  VARIATION  OF  THE  THRESHOLD  WITH  MCISTUR 
(COMPLETE  RESULTS) 
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conclusiows 

The  experiments  demonstrated  that  moisture  clearly  increases  the 
value  of  the  threshold  shear  velocity  of  sand  movement.  When  the  moisture 
has  its  origin  in  atmosphere  humidity,  the  variation  although  not  a  neg¬ 
ligible  one,  remains  rather  small  for  the  usual  range  of  air  humidities. 

On  the  contrary  when  the  water-content  of  the  sand  attains  the  values  of 
2  or  3%  the  wind  strength  necessary  to  initiate  the  movement  becomes 
considerable , 

If  w  is  the  water  content  expressed  in  percent,  the  threshold  shear 
velocity  U4  is  given  by, 

U#t  -  A  <1,8  +  0.6  log  «>  J  iZjl  gd 

When  atmostphere  is  responsible  for  the  sand  moisture,  it  is  preferable 
to  use  the  formula, 

0  a  A  ( 1  +  -  lL_ -  )  /  g  d 

*t  1  2  100  J  V  p  s 

where  h  is  the  air  humidity  expressed  in  percent.  In  both  formulae,  d 
represents  the  mean  grain  diameter,  o  and  p  are  the  density  of  sand  grains 
and  air,  respectively,  and  A  approximates  the  value  of  0,1  for  the  fluid 
threshold  and  0.08  for  the  impact  threshold. 
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